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The pericenter shift of a binary system represents a suitable observable to test for possible devia- 
tions from the Newtonian inverse-square law in favor of new weak interactions between macroscopic 
objects. We analyzed 13 years of tracking data of the LAGEOS satellites with GEODYN II software 
but with no models for general relativity. From the fit of LAGEOS II pericenter residuals we have 
been able to obtain a 99.8% agreement with the predictions of Einstein's theory. This result may 
be considered as a 99.8% measurement in the field of the Earth of the combination of the 7 and 
ft parameters of general relativity, and it may be used to constrain possible deviations from the 
inverse-square law in favor of new weak interactions parametrized by a Yukawa-like potential with 



strength a and range A. 
Earth radius. 
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Tests for Newtonian gravity and for a possible violation 
of the weak equivalence principle are strongly related and 
represent a powerful approach in order to validate Ein- 
stein's theory of general relativity (GR) with respect to 
alternative theories of gravity and to tune, from the ex- 
perimental point of view, gravity itself into the realm of 
quantum physics. Moreover, new long range interactions 
(NLRIs) may be thought of as the residual of a cosmolog- 
ical primordial scalar field related with the inflationary 
stage (dilaton scenario) Twenty-four years ago, the 
possibility of a fifth force of nature prompted new ex- 
perimental investigation of possible deviations from the 
gravitational inverse-square law Q. In fact, the devia- 
tions from the usual 1/r law for the gravitational poten- 
tial would lead to new weak interactions between macro- 
scopic objects. 

Interestingly, these supplementary interactions may be 
either consistent with Einstein's equivalence principle or 
not. In this second case, nonmetric phenomena will be 
produced with tiny, but significant, consequences in the 
gravitational experiments |3| • The feature of such inter- 
actions, which are predicted by several theories, is to pro- 
duce deviations for separations of masses ranging through 
several orders of magnitude, starting from the submil- 
limeter level up to the astronomical scale. Among the 
various techniques useful for the search of this additional 
physics to the various scales, the accurate measurement 
of the pericenter shift of binary systems may be used to 
test for a NLRI with a characteristic range comparable 
with the system semimajor axis 

These very weak NLRI are usually described by means 
of a Yukawa-like potential with strength a and range A 
and transmitted by a field of very small mass [i = h/Xc. 
If Goo represents the gravitational constant, and m s 



the mass of the primary body and of the satellite, r their 
separation, c the speed of light and h the reduced Planck 
constant, we can write: 
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where the strength a depends both on the mass-energy 
content of the sources and on their coupling strengths, 
Kq and K s , respectively. 

In the weak field and slow motion limit (WFSML) of 
GR, Einstein's equations reduce to a form quite similar to 
those of electromagnetism [5( . In Einstein's geometrody- 
namics and in the frame of a relativistic 3-body problem, 
where the two primaries are the Sun and the Earth and 
the test particle is represented by a satellite orbiting the 
Earth, the main precessions to which the satellite orbit, 
as a sort of enormous gyroscope, is subject to are com- 
monly known in the literature as: i) Einstein Q, ii) de 
Sitter (tJ, and iii) Lense-Thirring (LT) Q precessions. 
These precessions may be explained in terms of the ef- 
fects, on the orbital plane of the satellite, produced by the 
gravitoelectric and gravitomagnetic fields of the Earth, 
points i) and iii) , and by the effects arising from the cou- 
pling between the Earth's motion with the background 
field of the Sun, point ii). While Einstein's precession is a 
spin-independent secular effect, the other two precessions 
are usually interpreted as spin-orbit effects, in particu- 
lar, as frame-dragging effects, but with some differences: 
the de Sitter precession is frame-dependent, while the LT 
one is intrinsically related with the spin of the primary 
mass, i.e., with its rotation. Therefore, this precession 
must be related with intrinsic gravitomagnetism. We re- 
fer to Ciufolini and Wheeler and to Ciufolini [hJ for 
a deeper insight. 
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This Letter is devoted to showing some of our recent re- 
sults on the hrst simultaneous measurements of the cited 
rclativistic precessions in the field of the Earth using the 
two LAGEOS (LAser GEOdynamics Satellite) satellites. 
This work is new with respect to previous ones because 
we measure all the relativistic secular effects at one time. 
In particular, we focus on the satellites' pericenter secular 
advances 11], to which several non-Newtonian theories 
of gravity are sensitive. We analyzed 13 years of Satellite 
Laser Ranging (SLR) data of the two LAGEOS satellites 
using the NASA/GSFC software GEODYN II fl2|. This 
software is dedicated to satellite orbit determination and 
prediction, geodetic parameters estimation, tracking in- 
struments calibration, and many other applications in the 
field of space geodesy. The key ingredients of our mea- 
surement are (i) a consistent statement of the theory to 
be tested, (ii) the availability of a good test mass with re- 
lated high-quality tracking data and (iii) a modelization 
set for test mass dynamics and tracking. 

The first ingredient is far from being trivial: the rela- 
tivistic equations of motion can be formulated in princi- 
ple (due to general covariance) in whatsoever coordinate 
system; however, see Ashby and Bertotti (AB) [13(, a 
suitable choice of this system makes its physical interpre- 
tation clearer and simplifies its formulation. In their gen- 
eralized local inertial frame the main contribution to the 
test mass dynamics comes from the central body, while 
third-body effects show up only through tidal terms. The 
GR acceleration model included in GEODYN II follows 
the results of Huang et al. fl4j . and represents a gen- 
eralization of the AB model. The main feature of this 
model is that their noninertial geocentric frame retains 
all the merits of the inertial geocentric frame of AB, but 
it does not rotate with respect to the barycentric refer- 
ence frame. 

The second ingredient is given by the laser ranging 
data of LAGEOS satellites. The two are almost twins 
[H. LAGEOS, launched by NASA (1976), and LAGEOS 
II, launched by NASA/ASI (1992), have been designed 
spherical in shape, with high density and small area-to- 
mass ratio in order to minimize the effects of the subtle 
and complex nongravitational perturbations [l6j |. Their 
radius is just 30 cm and their mass about 407 kg fl7j . 
Their aluminum surface is covered with 426 cube-corner 
retro-reflectors for laser ranging from dedicated ground 
stations. The precision of the measurements is mainly 
related with the pulse width, which is usually rj 1 • 10~ 10 s 
down to 3T0 -11 s for the best laser ranging stations. The 
SLR data are available through the International Laser 
Ranging Service (ILRS) 18] in the form of normal points, 
with a root-mean-square (rms) down to a few mm, that 
corresponds to an accuracy in the orbit reconstruction at 
a few cm level, when using the best dynamical models. 
In our preliminary analyses we have been able to fit the 
orbit of the satellites at a 1-2 cm (rms) level in range. 

Regarding the third ingredient, the models included in 



GEODYN II are devoted to describe satellite dynamics, 
measurement procedure, and reference frame transforma- 
tions; they include [19]: (i) the geopotential (static and 
dynamic), (ii) lunisolar and planetary perturbations, (iii) 
solar radiation pressure and Earth's albedo, (iv) Rubin- 
cam and Yarkovsky-Schach effects (which need the satel- 
lite spin-axis coordinates), (v) SLR stations coordinates, 
(vi) ocean loading, (vii) Earth Orientation Parameters 
and (viii) measurement procedure. Usually, the models 
implemented in the code include the GR corrections in 
the parametrized post-Newtonian (PPN) formalism [2(| ■ 

In the analysis we performed, in order to solve for the 
relativistic secular precessions, we did not include in our 
setup such corrections. Moreover, we did not estimate 
any empirical accelerations as well as the satellites' ra- 
diation coefficient Cr, polar motion and universal time 
UT1 corrections, in such a way to avoid any possible ab- 
sorption of physical effects. Finally, in order to avoid the 
problems related with the spin modeling, also the Rubin- 
cam and Yarkovsky-Schach (YS) effects have not been in- 
cluded in our analysis. Concerning the estimated param- 
eters, besides the satellite state-vector, for each 15-day 
arc we estimated only measurement biases. The long-arc 
analysis of the orbits of geodetic satellites is a useful way 
to extract relevant information (i.e., model parameters) 
concerning the Earth's structure 21J . The physical infor- 
mation is concentrated in the satellite orbital residuals, 
that must be extracted from the orbital elements deter- 
mined during the fit (2^. The software fits the tracking 
data with all its models minimizing the difference be- 
tween the observed data and the computed ones, using a 
differential correction procedure. The final residuals are 
a measure of all the unmodeled effects, such as the GR 
ones, as well as of the poorly modeled and mismodeled 
effects and the noise in the tracking data. The relativis- 
tic precessions are effective in the orbital elements that 
define the orbit orientation in space, i.e., the orbit Eu- 
ler angles with respect to Earth's equatorial plane. These 
are the longitude Q of the orbit ascending node, the argu- 
ment of pericenter w, and the orbit inclination i. These 
elements are not equally sensitive, in their secular and 
periodic effects, to the relativistic precessions. Anyway, 
the argument of pericenter u) is sensitive to all of them, 
in particular to their secular effects. 

In our analysis we determined the orbital elements of 
the two LAGEOS satellites and then we computed the 
residuals with the method explained in Ref. [22j. This 
is the method developed by one of us in 1996, and that 
has been always used in the LT effect measurements per- 
formed so far [23] . This point represents a crucial aspect 
in this kind of analysis, because we need a reliable way 
to obtain the residuals in the orbital elements which re- 
tain the original concept of observed - computed quantity, 
which usually refers to the tracking observable, the range 
for SLR data. With regard to the background gravity 
field, in our setup we included two different models: (1) 
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EGM96 and (2) EIGEN-GRACE02S. The gravity field 
plays a very significant role in this kind of measurement. 
The uncertainties in its harmonic coefficients, especially 
in the even zonal ones, are the major source of system- 
atic effects, as we know very well in the case of the pre- 
vious measurements of the LT effect. The EGM96 model 
24( | , the current conventional model recommended by the 
International Earth Rotation Service (IERS) [l9j], is a 
multisatellite model derived over a time span of several 
years. The advantage of a multisatellite model resides in 
the different orbital characteristics of the satellites, such 
as different semimajor axis, eccentricity and inclination 
[25] . A clear disadvantage is represented by the fact that 
the tracking data and the quality of the orbit dynam- 
ical models are not homogeneous for all the satellites. 
EIGEN-GRACE02S [H has been derived by GRACE 
mission, and has the characteristic to improve the gravity 
field knowledge with a limited amount of data, in partic- 
ular in the medium and long wavelengths of its spectrum. 

In the following we focus on the results obtained from 
the analysis of the satellite's pericenter shift, in particu- 
lar, for LAGEOS II. Indeed, in the case of the pericen- 
ter the observable quantity is eAw; i.e., it depends on 
the satellite eccentricity e. Because LAGEOS II orbit 
is more eccentric, LAGEOS II is the best candidate for 
an accurate measurement of the total relativistic preces- 
sion of the pericenter. In Table Q] the results expected 
for the relativistic precession rates in the pericenter are 
shown. By inspection, the total relativistic precession of 
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FIG. 1. Integrated residuals (dots) of LAGEOS II pericenter 
with the linear trend (continuous line) fitted together with 
four periodic terms with the periods of the main spectral lines 
of the YS effect. The fractional discrepancy between the slope 
of the linear trend and the prediction of general relativity is 
just 2.8 ■ 1CT 4 . The starting epoch is MJD 49004. 



on the time span of the data analysis and on the num- 
ber of periodic effects which are fitted together with the 
linear term. The worst result we obtained, by changing 
some of the initial conditions and the number of adjusted 
parameters, has been a 0.21% discrepancy between the 
value of the recovered slope and the predicted one, see 
Table HI Therefore, for our analysis of the LAGEOS 



TABLE I. Rates in mas/yr of the secular relativistic preces- 
sion on the argument of pericenter of the two LAGEOS satel- 
lites (1 mas/yr = 1 milli-arc-second per year). 



Rates 


LAGEOS II 


LAGEOS 


Ad} E 


+ 3351.95 


+ 3278.77 


Auj lt 


- 57.00 


+ 32.00 


Au dS 


+ 10.69 


- 5.99 



the LAGEOS II pericenter is Aoj}f ~ 3305.64 mas/yr. 

The result of our 13 year analysis is shown in Fig- 
ure [1] for the satellite argument of pericenter advance 
in the case of the EIGEN-GRACE02S model. The Fast 
Fourier Transform (FFT) of the residuals in the peri- 
center rate confirms the presence of the unmodeled YS 
effect in the integrated residuals of Figure Q] We fit- 
ted the residuals with a linear trend plus four periodic 
terms [27]. These terms come from the FFT and cor- 
respond to the following main spectral lines of the YS 
effort [28]: A + lu (~ 257 days), X - u (~ 624 days), 
n + A + w (~ 485 days) and fi - A + u (~ 312 days) 
29] . For the linear trend slope we obtained a best value 
of Adjj} eas = 3306.58 mas/yr, that corresponds to a frac- 
tional discrepancy of about 0.03% with respect to the 
prediction of GR. The result of the fit mainly depends 



TABLE II. The fit results are mainly sensitive to the value of 
the intercept a of the fitting function. For a fixed to the value 
obtained from the residuals (280.9 mas), we obtained our best- 
fit result independently of the initial fixed value for the other 
parameters and the number of the periodic effects we added 
to the main four of the YS effect. Conversely, by varying the 
intercept up to ±50 mas, value that corresponds to the rms of 
the range residuals of the adjusted state vector, we obtained 
our worst result. Aa is the variation for the intercept while 
8a and 8b represent the adjustment of the given parameter. 



Aa [mas] 


8a [mas] 


Sb/ACj^f [%] 








± 0.03 


±10 


-9.9 


+ 0.06 


±20 


-19.9 


+ 0.10 


±30 


-29.9 


+ 0.13 


±40 


-39.9 


+ 0.16 


±50 


-49.9 


+ 0.21 



II pericenter general relativistic advance we assume the 
following conservative result: 



= 1 



- (0.28 ±2.14) • 10~ 3 . 



(2) 



The parameter e u may be considered at the post- 
Newtonian level, and measures possible deviations from 
Einstein's GR, where e w = 1. In the case of LAGEOS, we 
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obtained a worst best fit, with at least a 1% discrepancy 
with respect to the prediction of GR {Auj}" 1 ~ 3304.78 
mas/yr), due to the smaller eccentricity and the larger 
perturbations produced by the unmodeled YS effect [30J . 
The result obtained with the current analysis repre- 



produce a gravitomagnetic field, analogous to the mag- 
netic field p roduced by electric currents. 
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where also the contribution of the main systematic effects 
has been estimated. The secular effect is given by: 
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and it corresponds to the peak value at a range A = 
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consider our measurement as an upper bound for the 
strength of a possible long-range interaction; we obtain: 
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